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Abstract. We present H band (1.53/im-1.69/im) spectra of 57 isolated Be stars of spectral types 09-B9 and 
luminosity classes III,IV & V. The Hi Brackett (n-4) series is seen in emission from Br-11-18, and Fen emission 
is also apparent for a subset of those stars with Hi emission. No emission from species with a higher excitation 
temperature, such as Hell or Cm is seen, and no forbidden line emission is present. A subset of 12 stars show no 
evidence for emission from any species; these stars appear indistinguishable from normal B stars of a comparable 
spectral type. In general the line ratios constructed from the transitions in the range Br-11-18 do not fit case B 
recombination theory particularly well. Strong correlations between the line ratios with Br-7 and spectral type 
are found. These results most likely represent systematic variations in the temperature and ionization of the 
circumstellar disc with spectral type. Weak correlations between the line widths and projected rotational velocity 
of the stars are observed; however no systematic trend for increasing line width through the Brackett series is 
observed. 
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1. Introduction 

Be stars are defined as hot, non-supergiant stars that 
show, or have shown at some stage, emission lines in 
their spectrum. About 20% of B stars are Be stars hence 
they are an important part of the hot star population. 
Understanding them is crucial to our obtaining a complete 
picture of hot-star winds. They are rapid rotators (rotat- 
ing at a mean velocity of 70% of their break-up speeds, 
Porter 1996) and often have complex and variable emission 
line profiles that at sufficient resolution are double peaked. 
In addition they show an optical and infrared continuum 
excess. In spite of many attempts, no complete explana- 
tion for the Be phenomenon has yet been found. 

Observations have produced an empirical description 
of Be star circumstellar environments which is now gen- 
erally accepted: a dense "disc" exists in the equatorial 
plane (probably rotating in a Keplerian fashion), whilst 
over the polar regions there is a fast wind (velocities up 
to ^ 1000 — 2000km s""'^). The disc is ionized, and it is re- 
combination in the disc that gives the emission lines. The 
double peaked structure of these lines is then due to the 
velocity structure and to self-absorption in the disc. Gehrz 
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et al. (1974) showed that the optical and infrared excess 
in the systems could be explained by disc free-free emis- 
sion. Recent radio and optical interferometric data has 
confirmed the existence of the disc (Dougherty & Taylor 
1992, Quirrenbach et al. 1994, Stee et al. 1995). 

Studies of Be stars have traditionally concentrated 
on optical-near IR photometry and optical spectroscopy. 
However, the optical H i recombination lines do not func- 
tion as good diagnostics of the inner regions of the circum- 
stellar disc. Since it is likely that studying the inner re- 
gions of the disc will result in important constraints to the 
physical processes giving rise to the different wind regimes 
the need for observations of this region are pressing. Near- 
IR spectroscopy provides one tool to probe the innermost 
regions of the circumstellar disc. We have therefore ob- 
tained H and K band spectra of a sample of some ~60 Be 
stars, from B0-B9 to study this region of the circumstellar 
envelope for the first time. 

This paper is the third of a series on the optical and 
near IR spectral properties of a representative sample of 
58 Be stars. In Steele et al. (1999; Paper I) we discuss 
the the basic properties of our sample, such as spectral 
type, luminosity class and projected rotational velocity, 
and determine that no significant selection effects bias our 
sample. In Clark & Steele (2000; Paper II) we present 
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Fig. 1. H band spectra for Group 1 objects (I). 

the K band spectra of the sample and relate these to the 
underlying properties of the stars. In this paper we present 
H band (1.53-1.69 fim) spectra of 55 of the stars from 
Paper I, plus three additional objects (see Table 7). 

2. Observations & Data Reduction 

The sample of target objects contains objects from 09 to 
B8.5 and of luminosity classes III (giants) to V (dwarfs), 
as well as three shell stars. The sample was selected in 
an attempt to contain several objects that were typical of 
each spectral and luminosity class in the above range; it 
therefore does not reflect the spectral and luminosity class 
space distribution of Be stars, but only the average prop- 
erties of each subclass in temperature and luminosity. A 
spectral type and measure of v sin i was derived for each 
object in the sample and were presented in Paper I (only 
the spectral classifications are repeated here for sake of 
brevity) . The distribution of v sin i within each tempera- 
ture and luminosity class was carefully investigated and 
the conclusion drawn that there were no significant selec- 
tion effects biasing the average properties of the objects 
(see Paper I for details). 

The observations were carried out on the United 
Kingdom Infrared Telescope (UKIRT) on 1996 June 29 
and October 1-2 (see Tables 2-7), using the Cooled 
Grating Spectrometer (CGS4). The observations were 
made using the short focal length camera plus the 150 
line/mm grating, giving an ideal coverage from 1.53 
to 1.71/zm with a velocity resolution of ~ 100 km/s. 
Unfortunately, due to problems with the slit rotation 
mechanism, the wavelength coverage was slightly cur- 
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Fig. 2. H band spectra for Group 1 objects (II). 

tailed, leading to a final wavelength coverage of 1.53 to 

I. 69 /im, excluding the Hei 1.7004/im transition. Data re- 
duction was carried out in a similar manner as described 
in Paper II with correction for telluric features by ratioing 
observed G type standards with the solar spectrum to re- 
move features within their photosphere, and then ratioing 
the object spectrum with the modified standard spectrum. 

3. Discussion 

3.1. Line Identification 

The primary goal of these observations was to observe 
the Hi Brackett recombination series. The wavelength 
range chosen encompasses Br-18-11, as well as Fell lines 
at 1.534, 1.600 and 1.620/im (amongst others). Various 
He II transitions are also found within this wavelength re- 
gion, however given their absence in the K band spec- 
tra presented in Paper II, we do not expect to see these 
in emission. The spectra are presented in Figs. 1-6, and 
are presented in the same Groups we defined in Paper 

II, which were based on their K band spectral morphol- 
ogy (see Table 1). This approach has been adopted, rather 
than grouping the spectra on the basis of their optically 
derived spectral classifications since one aim of this work 
was to define a classification scheme based on near-IR 
spectroscopy alone, for use in the classification of heavily 
obscured stars, such as those often found in High Mass 
X-ray binary systems. 

Unsurprisingly, based on the features present in the K 
band spectra, the obvious lines seen in emission are Br- 
18-11, and Fe ii 1.576/um and 1.687/Ltm. A further feature 
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Table 1. K Band Spectral Groups and their corresponding optical spectral types (see Paper II for details) 



Group K Band appearance Spectral Type 



1 


Br 7 emission, He I features 


09e 


-B3e 


2 


Br 7 absorption, He I features 


09- 


-B3 


3 


Br 7 + Mg II emission, no He I features 


B2e 


- B4e 


4 


Br 7 absorption, no He I or Mg II features 


B4- 


B9 


5 


Br 7 emission, no He I or Mg II features 


B4e 


- B9e 



Table 2. Summary of measured Hi Brackett series, He I and Fell equivalent widths in A for Group 1 stars. Note 
that no correction is made in these tables for the underlying photospheric absorption line, and that a positive EW 
corresponds to emission, e indicates emission in the relevant transition that is to weak to measure (typically < 0.2A), 
the addition of a question mark represents an uncertainty as to its presence. A dash indicates that the feature was not 
within the wavelength range of that particular spectrum. The estimated error in the equivalent widths is 20 %. The 
spectral type of each star as determined in Paper I is given in column 2 (where an italicized type indicates that only 
a historical classification is available), and the observation date during 1996 is indicated in column 3. 



Object 


Spec. 


Obs. 


Brl8 


Brl7 


Brl6 


Brl5 


Fell 


Brl4 


Brl3 


Brl2 


Fell 


Brll 


Fell 




Type 


Date 










1.576 








1.679 




1.687 


CD -27 11872 


B0.5V-ni 


29/6 


8.3 


9.2 


9.8 


10.0 


0.8 


10.0 


9.7 


9.6 


e 


9.6 


e 


BD -13 893 


B1-3V 


1/10 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.2 


0.0 


1.4 




BD -12 5132 


BN0.2ra 


29/6 


2.8 


3.1 


3.5 


3.9 


1.0 


4.7 


4.6 


4.1 


e 


3.9 


e 


BD -1 3834 


B2IV 


29/6 


4.6 


5.3 


5.6 


6.2 


0.8 


6.3 


6.2 


5.7 


e 


5.6 


e 


BD +1 1005 


B1-3V 


2/10 


4.9 


5.0 


5.5 


6.3 


0.0 


6.1 


5.9 


5.6 


0.0 


5.7 




BD +4 1002 


B2-3III 


2/10 


0.2 


0.2 


0.3 


0.2 


0.0 


0.3 


0.7 


0.1 


0.0 


0.1 




BD +5 3704 


B2.5V 


29/6 


0.0 


0.0 


0.0 


0.1 


0.0 


0.2 


0.0 


0.0 


0.0 


-0.2 


0.0 


BD -1-29 4453 


B1.5V 


29/6 


8.7 


9.7 


9.9 


10.2 


0.9 


10.6 


11.5 


10.5 


e 


11.7 


e 


BD -1-36 3946 


BIV 


29/6 


4.6 


5.5 


5.5 


5.6 


0.4 


6.0 


5.7 


5.5 


e 


5.4 


e 


BD +45 3879 


B1.5V 


29/6 


7.8 


7.8 


7.8 


8.5 


0.7 


7.7 


8.4 


7.8 


e 


9.3 


e? 


BD +47 3985 


Bl-2sh 


29/6 


1.8 


2.1 


2.4 


1.9 


0.2 


1.8 


2.1 


2.4 


e 


4.1 


e 


BD +55 605 


BIV 


1/10 


0.7 


0.9 


1.4 


2.1 


0.0 


2.1 


2.4 


2.9 


0.0 


3.2 




BD +55 552 


B4V 


1/10 


0.1 


-0.4 


-0.1 


-1.2 


0.0 


-1.0 


-1.1 


-3.0 


0.0 


-0.3 




BD +56 469 


30-2111 


1/10 


-0.4 


-0.4 


-0.3 


-0.3 


0.0 


-0.3 


-0.3 


-0.6 


0.0 


-0.6 




BD +56 473 


BIV-III 


2/10 


3.4 


3.5 


3.7 


3.9 


0.5 


4.0 


3.9 


3.8 


e 


4.5 




BD +56 478 


B1.5V 


1/10 


0.9 


1.2 


1.4 


1.5 


0.0 


1.5 


1.4 


1.9 


0.0 


2.1 




BD +56 511 


BlIII 


1/10 


0.5 


0.4 


0.7 


1.0 


0.0 


0.1 


0.3 


0.7 


0.0 


0.3 




BD +56 573 


B1.5V 


1/10 


5.9 


7.2 


7.8 


8.6 


2.1 


9.8 


8.8 


9.7 


e 


10.3 




BD +57 681 


B0.5V 


1/10 


0.6 


0.4 


0.9 


1.2 


0.0 


0.9 


1.0 


1.6 


0.0 


2.0 




BD +58 2320 


B2V 


29/6 


2.9 


3.6 


4.1 


4.0 


0.0 


3.2 


4.6 


3.8 


0.0 


5.9 


0.0 



Table 3. Summary of H I Brackett series for Group 2 stars. See Table 2 for explanation 



Object 


Spec. 


Obs. 


Brl8 


Brl7 


Brl6 


Brl5 


Fell 


Brl4 


Brl3 


Brl2 


Fell 


Brll 


Fell 




Type 


Date 










1.576 








1.679 




1.687 


CD -25 12642 


B0.7ni 


29/6 


-0.3 


-0.5 


-0.4 


-0.8 


0.0 


-1.7 


-2.9 


-3.0 


0.0 


-4.0 


0.0 


BD +20 4449 


BOIII 


29/6 


0.0 


0.0 


0.0 


-0.5 


0.0 


-1.0 


-1.7 


-2.0 


0.0 


-3.0 


0.0 


BD +25 4083 


B0.7in 


29/6 


-0.4 


-0.2 


-0.1 


-0.7 


0.0 


-1.5 


-2.1 


-2.5 


0.0 


-3.2 


0.0 


BD +28 3598 


09n 


29/6 


0.0 


0.0 


0.0 


0.0 


0.0 


-0.3 


0.3 


-0.7 


0.0 


-1.5 


0.0 


BD +29 3842 


Bin 


29/6 


-0.3 


-0.3 


-0.3 


-0.8 


0.0 


-1.4 


-1.6 


-2.3 


0.0 


-2.3 


0.0 


BD +37 3856 


B0.5V 


29/6 


0.0 


0.0 


0.0 


-0.5 


0.0 


-0.8 


-1.7 


-2.3 


0.0 


-2.8 


0.0 


BD +45 933 


B1.5V 


1/10 


-0.5 


-0.9 


-1.1 


-1.6 


0.0 


-1.8 


-1.3 


-2.0 


0.0 


-1.9 




BD +56 493 


BIV 


1/10 


-0.3 


-0.1 


-0.2 


-0.7 


0.0 


-1.3 


-2.3 


-1.7 


0.0 


-2.9 





in the blue wing of Br-11 (A '^1.678/im) is also seen in a 
subset of spectra, the identity of which is discussed further 
below. No emission from species of higher excitation (such 
as Hen, Niii or Cm) was observed. The presence and 
equivalent widths (EW) of these features are summarised 



in Tables 2 7. Equivalent widths and FWHM were mea- 
sured using the ABLINE routine of the FIGARO software 
package. This measures EW by integrating the line flux 
relative to an interpolated continuum produced by poly- 
nomial fitting in the vicinity of the line - it therefore makes 
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Table 4. Summary of H I Brackett series for Group 3 stars. See Table 2 for explanation 



Object 


Spec. 


Obs. 


Brl8 


Brl7 


Brl6 


Brl5 


Fell 


Brl4 


Brl3 


Brl2 


Fell 


Brll 


Fell 




Type 


Date 










1.576 








1.679 




1.687 


BD -8 929 


B2V 


2/10 


4.6 


4.5 


4.5 


4.6 


0.0 


5.0 


5.2 


4.6 


0.0 


4.9 




BD +42 4538 


B2.5 


29/6 


8.2 


8.7 


9.2 


8.4 


0.0 


9.0 


9.3 


9.2 


0.0 


10.6 


0.0 


BD +47 183 


B2.5V 


29/6 


5.6 


5.6 


5.7 


5.8 


0.3 


5.6 


6.4 


6.3 


e 


7.0 


e 


BD +47 857 


B4IV 


1/10 


2.4 


3.5 


4.3 


4.4 


1.5 


4.4 


4.5 


5.3 


e 


5.7 




BD +47 939 


B2.5V 


1/10 


3.6 


4.1 


4.3 


4.6 


0.9 


5.4 


5.2 


5.2 


e 


6.4 





Table 5. Summary of H i Brackett scries for Group 4 stars. See Table 2 for explanation 



Object 


Spec. 


Obs. 


Brl8 


Brl7 


Brl6 


Brl5 


Fell 


Brl4 


Brl3 


Brl2 


Fell 


Brll 


Fen 




Type 


Date 










1.576 








1.679 




1.687 


BD -19 5036 


B4III 


29/6 


-0.8 


-1.4 


-2.0 


-2.5 


0.0 


-4.8 


-5.4 


-6.3 


0.0 


-7.5 


0.0 


BD +17 4087 


B6III-V 


29/6 


-1.2 


-1.4 


-2.1 


-2.9 


0.0 


-3.5 


-5.1 


-5.2 


0.0 


-5.7 


0.0 


BD +19 578 


B8V 


1/10 


-0.3 


-0.4 


-0.4 


-1.2 


0.0 


-2.1 


-4.1 


-5.3 


0.0 


-5.9 




BD +30 3227 


B4V 


29/6 


0.0 


0.0 


0.0 


-1.0 


0.0 


-3.4 


-5 


-5.7 


0.0 


-7.5 


0.0 



Table 6. Summary of H I Brackett series for Group 5 stars. See Table 2 for explanation. 



Object 


Spec. 


Obs. 


Brl8 


Brl7 


Brl6 


Brl5 


Fell 


Brl4 


Brl3 


Brl2 


Fen 


Brll 


Fell 




Type 


Date 










1.576 








1.679 




1.687 


CD -27 13183 


B7V 


29/6 


0.0 


0.0 


-1.0 


-1.4 


0.0 


-2.0 


-3.1 


-3.2 


0.0 


-4.80 


0.0 


BD -20 5381 


B5V 


29/6 


1.0 


1.0 


1.7 


1.1 


0.0 


0.7 


1.4 


0.7 


0.0 


0.8 




BD -0 3543 


B7V 


29/6 


0.0 


0.0 


0.0 


-1.2 


0.0 


-2.2 


-2.4 


-2.7 


0.0 


-3.3 


0.0 


BD +0 1203 


B5III 


29/6 


1.1 


0.7 


0.7 


0.8 


0.0 


0.9 


0.8 


1.2 


0.0 


2.0 


0.0 


BD +2 3815 


B7-8sli 


29/6 


-0.5 


-0.8 


-1.2 


-2.0 


0.0 


-3.7 


-3.5 


-4.1 


0.0 


-4.6 


0.0 


BD +21 4695 


B6III-V 


29/6 


0.0 


-0.3 


-0.8 


-1.3 


0.0 


-2.8 


-3.4 


-5.4 


0.0 


-6.4 


0.0 


BD +27 3411 


B8V 


29/6 


0.0 


-0.4 


-0.7 


-1.1 


0.0 


-2.4 


-4.5 


-5.6 


0.0 


-6.5 


0.0 


BD +37 675 


B7V 


1/10 


1.2 


0.6 


1.0 


0.0 


e 


-0.2 


-1.7 


-2.0 


e 


-1.5 




BD +43 1048 


B6IIIsli 


1/10 


0.2 


0.2 


0.2 


0.3 


0.0 


0.5 


0.7 


1.4 


0.0 


1.2 




BD +46 275 


B5III 


2/10 


-0.1 


-0.9 


-2.2 


-2.4 


0.0 


-4.9 


-4.3 


-5.2 


0.0 


-6.0 




BD +49 614 


B5ni 


1/10 


0.2 


0.3 


0.6 


0.5 


0.0 


-1.6 


-3.3 


-2.7 


0.0 


-2.7 




BD +50 825 


B7V 


2/10 


0.3 


0.3 


0.2 


-1.0 


0.0 


-0.2 


-1.7 


-2.0 


0.0 


-1.5 




BD +50 3430 


B8V 


29/6 


0.5 


0.5 


0.6 


-1.7 


0.0 


-1.7 


-1.9 


-2.6 


0.0 


-4.5 


0.0 


BD +51 3091 


B7III 


29/6 


-1.0 


-1.8 


-2.4 


-3.8 


0.0 


-4.3 


-6.2 


-6.8 


0.0 


-6.7 


0.0 


BD +53 2599 


B8V 


29/6 


-1.1 


-1.9 


-1.9 


-2.7 


0.0 


-2.4 


-3.7 


-5.1 


0.0 


-7.5 


0.0 


BD +55 2411 


B8.5V 


29/6 


-0.4 


-0.6 


-1.2 


-2.1 


0.0 


-4.0 


-5.0 


-6.0 


0.0 


-7.5 


0.0 


BD +58 554 


B7V 


1/10 


0.0 


0.8 


1.2 


0.7 


0.0 


0.2 


0.7 


0.6 


0.0 


0.3 





no assumption about line shape. FWHM are measured by 

looking at the line width which encompasses 68% of the 
line flux and multiplying by 1.18. For a line of Gaussian 
profile this gives exactly the FWHM, and for any reason- 
ably centrally concentrated profile (as ours are) it gives a 
result very close to the FWHM that is more robust than 
simple Gaussian fitting. 

Brackett scries emission is seen in all stars of Groups 
1, 3 and 5 (which all possess Br7 emission), with the pos- 
sible exception of BD +51 3091, a B7 HI star (Fig. 6). Of 
Groups 2 and 4, which were defined on the basis of the 
absence of Br7 emission in the K band spectra, only BD 
+30 3227 appears to show evidence for partial infilling of 
Br- 11 14 (Fig. 3), while the remainder appear to possess 
pure absorption spectra. 



Fe II 1.576/Lim emission is seen in a total of 13 stars, 
compared to a total of 19 stars showing Fell 2.089/im 
emission. Of the 33 stars where the wavelength coverage 
encompassed Fell 1.688/Lim, the line was in emission in 6 
stars, 5 of which are Group 1 objects, the final star, MWC 
659, being one of the 3 objects without corresponding K 
band observations (see Tables 2-7). We note that all stars 
with Fe II 1.576/im emission also show Fc ii 1.688/im emis- 
sion. The correlation between the presence of Fe ii emis- 
sion in the H and K bands is weaker, with a total of 8 stars 
showing Fc II emission in either the H or K band, but not 
in the other. However, as with the K band spectra, Fe ii 
emission only occurs in the H band spectra of those stars 
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Table 7. Summary of H I Brackett series for new stars. See Table 2 for explanation. 



Object 


Spec. 


Obs. 


Brl8 


Brl7 


Brl6 


Brl5 


Fell 


Brl4 


Brl3 


Brl2 


Fen 


Brll 


Fell 




Type 


Date 










1.576 








1.679 




1.687 


BD +54 2348 


B2V 


29/6 


6.2 


7.0 


7.3 


7.9 





7.6 


7.7 


7.9 


0.0 


8.7 


0.0 


BD +54 2718 


B2 III 


29/6 


-0.5 


-0.4 


0.0 


0.0 


-0.7 


-1.0 


-0.9 


-1.4 


0.0 


-3.0 


0.0 


MWC 659 


BOIIIpe 


29/6 


7.8 


8.9 


9.8 


9.4 


0.8 


10.7 


10.3 


10.9 


e 


11.3 


e 



BD +55 605 




1.55 



1.6 

Wavelength (microns) 



1.65 



Group 4 



1.55 



1.6 

Wavelength (microns) 



BD +37 3856 
BD +56 493 
BD +45 933 
BD +20 4449 
BD +25 4083 
CD -25 12642 
BD +29 3842 
BD +28 3598 



1.65 



Fig. 3. H band spectra for Group 1 (top panel) objects 
(III), Group 2 (middle panel) objects and Group 3 objects 
(bottom panel). 



with strong Brackett series emission; 

(EWBr7 >8A). 

Two possible identifications exist for the emission fea- 
ture at ~1.678Atm; Fell 1.679Atm and [Fell] 1.678/im. Of 
these, we favour Fell 1.679/xm since the feature is only 
seen in stars that show emission in the other two Fc ii 
lines, and no emission is seen in the (albeit weaker) [Fe ii] 
transitions at 1.534/im and 1.644/zm (Hamann & Persson 
1989). 

Overall, as regards spectral classification using the H 
band spectra of Be stars, it is apparent that with the ab- 
sence of many pure photosphcric features that arc uncon- 
taminated by disc emission it is only possible to perform 
a very general spectral classification into "early" (BOe- 
B4e) and "late" (B5c B9c) spectral types. The problem 
of the lack of photosphcric features uncontaminated by 
emission from the circumstellar disc is compounded by 
the (('xp(H tccl) lack of emission from species with a wide 



^ Note that in this paper we will employ the convention that 
positive equivalent widths indicate emission features. 



Fig. 4. H band spectra for Group 4 objects. 



range of excitation energies (although we speculate that 
the presence of He 1 1.700/im emission may function as an 
additional, valuable diagnostic of early spectral types). 

3.2. Brackett Line Strengths 

In order to make an accurate comparison of Brackett line 
strengths between objects of different spectral type, it is 
necessary to remove the effect of the underlying photo- 
spheric absorption lines from the spectra. To carry out 
this correction we used the equivalent widths for normal 
(non-emission) B stars presented for Br-7 by Hanson et al. 
(1996) and for Br 11-18 by Steele & Clark(2001). The cor- 
rection itself simply consisted of subtracting the equivalent 
width derived from a least squares fit to the appropriate 
transition for the normal B stars from the emission line 
equivalent width. The discussion and figures in this section 
therefore refer to this corrected equivalent width. We note 
that within the small wavelength range encompassed by 
Br- 11 to Br-18, these corrected equivalent widths may be 
divided by one another and treated as flux ratios without 
introducing significant error. 
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Fig. 5. H band spectra for Group 5 objects (I). 
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Fig. 6. H band spectra for Group 5 (upper panel) objects 
(II) and new objects (lower panel). 



Wc plot the EW of Br-11, 14 & 18 against spectral type 
in Fig. 7 (the remainder of the lines are not plotted due 
to reasons of space). As with Br7 (Paper II) we find that 
while a linear correlation between spectral type and EW 
is absent, a strong trend in the upper envelope of the line 
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Fig. 7. Equivalent width (EW) of Br-11 (top panel), Br- 
15 (middle panel ) and Br-18 (bottom panel) in A (where 
positive EW's indicate emission) against spectral type. 
The equivalent widths have been corrected for the un- 
derlying photospheric absorption as indicated in the text. 
Triangular symbols represent luminosity class III, squares 
luminosity class IV and circles luminosity class V. 

strengths is present, with lower mean and maximum line 
strengths for the later (~B5-B9) spectral types. We find 
no evidence of systematic diff'erences between the emission 
characteristics of stars of differing luminosity classes. 

In Fig. 8 we plot the EW of Br-11 against those of 
Br-12, -15 and -18, while in Fig. 9 we plot the EW of Br- 
18 versus Br-12, -15, -16 and -17. Note that the empty 
symbols cluster around EW of zero in both papers, with 
a typical deviation of ~ 1 A. As these are non-emission 
line objects, this indicates that the accuracy of our correc- 
tion for the photospheric absorption component is of that 
order. We also note that there appears to be no signifi- 
cant difference in these graphs between the various spec- 
tral groups. 

Given that the Br 11-18 lines arise from high levels in 

the hydrogen atom they might be expected to be optically 
thin, and therefore well represented by case B recombina- 
tion theory. In order to test this, we plot as a straight 
line in Fig.s 8 and 9 the expected case B line ratios from 
Storey & Hummer (1995) for ionized hydrogen. The pre- 
dicted line ratios are constant to within a few per-cent for 
temperatures in the range 7500 - 40000 K and electron 
densities in the range 10^° — 10^'* cm~^, and are therefore 
should not be sensitive to these conditions within the disk. 
From the figures it is apparent that although there is rea- 
sonable agreement between case B theory and observation 
in the line ratios between closely adjacent lines (e.g. Br-11 
and Br-12, or Br- 17 and Br-18), for more separated tran- 
sitions the case B ratios do not seem to fit the data well. 
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Fig. 8. EW of Br-11 against Br-12 (top panel), Br-15 
(middle panel) and Br-18 (bottom panel). Filled circles 
represent objects from Groups 1, circles Group 2, stars 
Group 3, empty squares Group 4 and filled squares Group 
5 (see Table 1 for group definitions). The line indicates 
the expected line ratio for case B recombination. 
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EW(Br18) 

Fig. 9. EW of Br-18 against Br-12, -15, -16 and -17. 
Symbols and line as Fig. 8. 

This may reflect either a complex disk temperature and 
density profile (as discussed in the following paragraph) or 
that the optically thin and/or LTE assumptions for these 
lines are not valid. Spectra extending further down the se- 
ries (e.g. to Br-22) where the lines are weaker would help 
to resolve this question. 



Fig. 10. Plot of the EW of Br7 (data from Paper II) 
against Br-11 (top panel), Br-15 (middle panel) and Br- 
18 (bottom panel; aU EW in -A). Symbols as Fig. 8. 

In Fig. 10 we plot the EW of Br7 (taken from Paper 

II) against that of Br-11, 15 & 18. Unlike the higher tran- 
sition ratios of Fig. 2 and 3 here there arc large differ- 
ences in the flux ratios between the groups, with a pro- 
gression in the ratio Br-15/Br7 from ~ 0.5 (Group 1), 
through ~ 0.2 (group 3) to ~ 0.0 (Group 5). This effect 
is a real reflection of differences in line flux and not just 
due to continuum differences over the larger wavelength 
range for the Br-15/Br7 equivalent width ratio. At first 
sight it may seem that the stronger circumstellar free-free 
excesses of the Group 1 objects could account for this ef- 
fect, by reducing the Br7 EW for these objects. However 
this effect is small. The two continuum regions are well 
approximated by the photometric H and K bands, and 
Howells et al. (2001) show that the mean E{H — K)cs 
only varies from ~ 0.12 for Group 1 objects to ~ 0.03 
for Group 5 objects, insufficient to cause the observed 
EW ratio changes. We also note that this effect is some- 
what mitigated by the change in intrinsic colours between 
Groups 1 {{H - K)o - -0.04) and 5 {{H - K),, - -0.01) 
(Koornneef 1983). We therefore believe that this progres- 
sion is likely to be due to systematic changes in the tem- 
perature, degree of ionization, and structure within the 
discs as the temperature and flux of the underlying stars 
changes. Marlborough et al. (1997) present simulations of 
selected near-IR Hi transitions for both isothermal discs, 
and discs with a simple radial temperature gradient, and 
find that varying the disc temperature does indeed lead 
to changes in the line decrement. However, recent work 
by Millar & Marlborough (1998) shows that these simple 
models for the disc temperature are incorrect, and that 
the equilibrium temperatures for Be star discs are a com- 
plex function of disc radius, particularly in the inner disc 
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Table 8. Summary of best fits between the FWHM of 
Brll-18 and the mini of the individual star. The best fit 
from Br7, taken from Paper II is also shown for compari- 
son. 
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Fig. 11. Plot of the FWHM (in A) of Br-11, 12, 15 and 18 
plotted against projected rotational velocity. Dotted line 
indicates FWHM=tsini from Paper I. 



regions responsible for near-IR recombination line emis- 
sion. Further evidence for this is found in our inability to 
fit case B line ratios to the higher transitions as noted 
previously. Consequently, we defer detailed discussion of 
the systematic variations in the Brackett line decrement 
for a future paper, where modeling of the line fluxes will 
be accomplished for the full optical-IR dataset. 

3.3. Line widths and projected rotational velocities 

In Paper II we found correlations between the full width 

half maximum (FWHM) of Br7 and He I 2.058/im lines 
and the projected rotational velocity of the stars, in a sim- 
ilar fashion that those found for Ha, 0, and selected opti- 
cal Fc II lines by Hanuschik (1989). He further foimd that 
the intrinsically weaker lines such as the higher Balmer 
lines and Fe ii were systematically broader than Ha which 
he attributed to the weaker lines being produced in a 
smaller region of a quasi- Keplerian disc than the in- 
trinsically stronger Ha line, where the rotational velocity 
would be greater, and thus the line broader. Simulations 
by Marlborough et al. (1997) also demonstrate the same 
effect, as do the longer wavelength infrared spectra of 7- 
Cas presented by Hony et al. (2000). 

In Fig. 11 we plot the full width half maximum 
(FWHM) of Br-11, 12, 15 and 18 for stars of Groups 1 
and 3 against their projected rotational velocities from 
Paper 1 (stars from Group 5 were excluded due to the 
bias introduced into the FWHM by the underlying pho- 
tosphcric feature). Weak correlations were found between 
the FWHM and the projected rotational velocities by ap- 
plying Spearman's Rank Correlation to the datasets, and 
the results of least square fits to the data arc presented in 
Table 7. It is interesting to note that the minimum mea- 
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sured FWHM is aroimd 300 km/s (considerably larger 
than the instrumental resolution of 100 km/s). This 
means that for the relatively slow (vsini ~ 100 — 200 
km/s) rotators, the FWHM for these lines is often more 
than twice the v sin i of the underlying star, possibly in- 
dicating more rapid rotation at the inner edge of the disk 
than at the surface of the star. We also note three Group 
1 (BD +47 3985, BD -^58 2320 and BD +1 1005) and 
one Group 3 (BD -|-42 4538) appear to have very large 
Br-11 FWHM of greater than around 500 km/s. However 
these objects also show double peaks in their line profiles, 
making the derivation (and interpretation) of FWHM un- 
certain. 

In Fig. 12 we plot the FWHM of Br-18, 15 & 12 vs 
Br-11. It can be seen that there is no trend to broader 

line profiles with higher Brackett series transitions as was 
suggested by the results of Hony et al. (2000). We attribute 
this result primarily to the narrow wavelength (and hence 
transition) range of our spectra compared to the work of 
Hony et al. (2000) whose spectra extend to the limit of 
each series. 



4. Conclusions 

We have analysed the H band (1.53/im 1.69/im) spectra of 
61 Be stars and found emission from H i Brll-18 and Fe 11 
1.576 & 1.698/im. A further emission feature at 1.679/xm 
is present in a subset of 14 stars; this is also likely to be 
Fell emission. The Br 11-18 line ratios of non adjacent 
lines are not well fit by the case B assumption, and strong 
systematic trends in both line strength and the ratio of 
the higher Brackett series strengths to Br-7 with spectral 
type (but not luminosity class) are observed. This is likely 
to be due to systematic changes in the temperature and 
degree of ionization between the discs of stars of early (BO- 
B4) and late (B5-B9) spectral type. We find that analysis 
of H band spectra alone only allow the classification of 
stars into "early" (B0e-B4e) or late (B5e-B9e) types; no 
determination of the luminosity class of the object can 
be made. This is due to the lack of any uncontaminatcd 
photospheric features in this region, the lack of emission 
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Fig. 12. FWHM of Br-18, 15 & 12 vs Br-11 for the ob- 
jects showing emission features. No systematic difference 
between the FWHM of Br-11, and those of the higher 
transitions is apparent. The larger scatter of points in the 
lower panel is likely to be due to the weakness of the Br-18 
line, and the consequent difficulties in obtaining an accu- 
rate measurement of the FWHM. Symbols as Fig. 8. 

features encompassing a wide range of excitation temper- 
atures and the similarity of Brackett line strengths and 
decrements for the early type stars. As with Br7, we find 
weak correlations between the FWHM of Brll-18 and the 
projected rotational velocity of the midcrlying stars. We 
find no systematic trend in FWHM through the Brackett 
series down to Br-18. 
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